Orbital magnetic response and the anisotropy of magnetic susceptibility in the 

Iron-based superconductors 
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We propose that the orbital angular momentum of the conduction electrons in the Iron-based 
superconductors is activated in their low energy physics. Using a five-band tight-binding model 
derived from fitting the LDA band structure, we find that the orbital magnetic susceptibility of 
the conduction electrons in such a multi-orbital system is several times larger than the Pauli spin 
susceptibility and is comparable in magnitude to the observed total magnetic susceptibility. The 
orbital magnetic susceptibility in the Fe-As plane(xi,) is found to be larger than that perpendicular 
to the Fe-As plane(xl,) by a factor about two and the total magnetic susceptibihty in the normal 
state can be fitted with formula x{T, d) ~ Xs{T) + XL{d), where Xa{T) is the temperature dependent 
isotropic part due to spin and xl{0) is the temperature independent anisotropic part due to orbital. 
In the superconducting state, Xl is found to be significantly reduced as the pairing gap develops, 
while Xl is almost not afi'ected by the superconducting transition. We argue the large anisotropy 
observed in the bulk magnetic susceptibility and the Knight shift in the Iron-based superconductors 
should be attributed to the orbital magnetic response of their conduction electrons. 

PACS numbers: 



One special feature of the newly discovered Iron-based 
superconductors^"— is their multi-orbital nature. In con- 
ventional superconductors, in which only one orbital 
plays an essential role around the Fermi energy, the or- 
bital angular momentum is quenched either as a result 
of the s-wave character of the orbital, or by the crystal 
field splitting effect. However, from LDA band struc- 
ture calculation, it is reported that all the five Fe 3d 
orbital play essential role in forming the low energy de- 
gree of freedom around the Fermi energy and the crys- 
tal field splitting is extremely small. This situation has 
caused great complexities in the model study. However, 
it also generates the interesting opportunity to explore 
the physics of the orbital degree of freedom in this sys- 
tem. The orbital character near the Fermi energy has 
now been extensively studied by angle-resolved photone- 
mission spectroscopy (ARPES) and the LDA result is to 
a large extent confirmedi^"- 

Many novel properties of the Iron-based superconduc- 
tors have been attributed to their multi-orbital nature. 
For example, the intimate relation between the structural 
and magnetic phase transition o^'^'^ has been proposed to 
originate from an orbital-related mechanismjiiii^. The 
same picture also provides a reasonable interpretation 
for the unusual in-plane anisotropy in resistivit y "'^^'^^ and 



tributed to the response of its spin degree of freedom, 
which in the absence of the spin-orbital coupling is 
isotropic. The orbital magnetic susceptibility, which is 
in principle anisotropic, is usually small as a result of 
the quenching of the orbital angular momentum of the 
conduction electrons. The small remnant orbital mag- 
netic response, namely the well known Van Velck para- 
magnetism, is controlled by the large gap separating the 
occupied bands from the unoccupied atomic levels that 
carry orbital angular momentum and is expected to show 
negligible temperature dependence. However, if the con- 
duction electron itself carries orbital angular momentum, 
the orbital magnetic response of the system would be 
much larger and be sensitive to the changes of electronic 
state of the system. 

Measurement of the magnetic susceptibility on the 
Iron-based superconductor has produced several intrigu- 
ing results. Firstly, the magnetic susceptibility is found 
to show significant temperature dependence J^"— Such 
a temperature dependence is unexpected for a weakly 
correlated system with a large band width. In the lit- 
erature, this unusual temperature dependence is either 
attributed to the strong correlation effect, or to the prox- 
imity of the system to a semimetal phase^^r^. Sec- 
ondly, the magnetic susceptibility is found to be strongly 
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the dxz and dyz band splitting observed in the ARPES anisotropic.—"— As shown in Fig. [l] the susceptibility in 

the Fe-As plane is much larger than that perpendicular 
to the Fe-As plane and their difference is nearly temper- 
ature independent. An understanding of such anisotropy 
in the magnetic susceptibility is still absent and we will 
show that it can be originated from the contribution of 
the orbital angular momentum of the conduction elec- 
trons. 



measurementjS In these situations, the multi-orbital na- 
ture manifests itself in the form of a static structure, a 
more interesting possibility is that the orbital degree of 
freedom appearing as a dynamic mode in the low energy 
physics, contributing to various kinds of response and re- 
laxation processes. The purpose of this paper is to inves- 
tigate the contribution of the orbital angular momentum 
of the conduction electrons to the magnetic susceptibility. 

The magnetic susceptibility of a metal is usually at- 



In this paper, we study the orbital magnetic response 
of the conduction electrons in the Iron-based supercon- 
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FIG. 1: The experimental results of the magnetic susceptibil- 
ity(a) and Knight shift(b) for some Iron-based supercondu- 
tors. The magnetic susceptibility data of Bal22 system are 
quoted from Ref. ,18, data Cal22 system from Ref. [iTlfboth 
in unit of 10~^ emu/Oe), and data of Srl22 system from 
Ref. |li(in unit of 10"^ emu/mol). The Knight shift data 
are quoted from Ref. [23. Here the superscripts xy/x and z 
denote the direction of the applied magnetic field. 



ductors. It is found that the orbital magnetic susceptibil- 
ity of the conduction electrons is several times larger than 
the Pauli spin susceptibility and is comparable in magni- 
tude with the observed total magnetic susceptibility. The 
orbital magnetic susceptibility in the Fe-As plane(xf,) is 
found to be larger than the susceptibility perpendicular 
to the Fe-As plane(x|^) by a factor about 2. In the nor- 
mal state, the magnetic susceptibility can be separated 
into a temperature dependent isotropic part Xs {T) due to 
spin response and a temperature independent anisotropic 
part Xi(^) due to orbital response, where is angle be- 
tween the magnetic field and the normal of the Fe-As 
plane. In the superconducting state, when a s^-wave 
intraband pairing is assumed, Xl found to be signifi- 
cantly reduced as the pairing gap develops, while Xl 
almost unaff'ected. Unlike the Pauli spin susceptibility, 
the orbital magnetic susceptibility is nonzero in the zero 
temperature limit as a result of the interband contribu- 
tion to the orbital magnetic response. 

The properties of the Iron-based superconductors are 
sensitive to model parameters. To understand its physics, 
a realistic band structure is indispensable. In this paper, 
we will adopt the five-band tight-binding model derived 
from fitting the LDA band structurei^. The model reads 
(following the notations of Ref. [isl ). 



3(i3z2_r2, 3dxz, 'idyz, 'idxy and 'id^^-yi respectively. 
t(Aa;, Ay; /i, v) denotes the in-plane hopping integral be- 
tween the /j,-th and v-ih orbitals at the lattice distance 
(Aa;, Ay), is the on-site energy of the /x-th orbital. The 
values of these model parameters are given in ReflT5l. 
The model Hamiltonian Eq.([T]) can be diagonalized in 
momentum space, in which it takes the form 



k,Q:,cr 



(2) 



where ek,a is the band energy for the a-th band. 

To describe the superconducting state, we model 
the pairing potential with a sign-changing s-wave form 
AQ,^(k, T) = A(T)5c</3 cos(fc3;) cos(fcj,). Note that only 
the intraband pairing is considered. Such a pairing 
symmetry is consistent with the weak coupling spin- 
fluctuation-exchange mechanism^^ or the strong coupling 
superexchange mechanism^*' for the superconductivity in 
the Iron-based superconductors. However, as will be 
clear below, the detailed form of the pairing potential 
(except its intraband pairing nature) is not essential for 
the uniform susceptibility. So we will take the pairing 
potential adopted only as a simplified way to induce a 
full gap on the Fermi surfaces of the system. 

The orbital magnetic susceptibility is defined through 
the correlation function of the orbital magnetic moment 
in the following way 

X2(q,r) = -(5LMB)'(T.L'^(q,r)L''(-q,0)), (3) 

in which L°(q, r) denotes the Fourier component of the 
orbital magnetic moment density in the a-direction. The 
orbital magnetic moment on a given site i is defined as 
Li = T.yy,s4,v,sKy<^^,'^',s, where l^y is the matrix 
element of the orbital magnetic moment in the space 
spanned by the five orbital 3c?322_r2, 'idxz^ idyz^ idxy 
and 3(^2,2 _y2. Here we set gLfJ'B = 1- 

To derive an expression for the matrix element , 
we approximate the Wannier functions with the Fe 3d 
atomic orbital. We then have 
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where \fi) {ji = 1, 2, 3, 4, 5) denote the five Wannier func- 
tions, \l,m) oc Y;™ denotes the eigenstate of the orbital 
angular momentum (P, Z^), whose eigenvalue for P is m. 
With the transformation Eq.(|3]), the matrix ^, and j,, 
can be found as 



where fijV = 1,2,3,4,5 denote the five Fe 3d orbital 
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The orbital magnetic susceptibility is readily obtained 
as follows 



XUT) = 



1™ TT 



1 



^'k.Q — £^k+q,^ 
"^k.Q^k+q,^ + Ak,QAk+q,/3 



(5) 



£^k,a-Ek+q,/3 
1 - /(^^k+q,^) - fiEk^a) 



1 



Ek,a + -Ek+q,/3 
Ck,QCk+q,/3 + Ak,aAk+q,/3 



\^W,a,l3\ 



Pk,Q,/3| 



£^k,Q£^k+q,,9 

Here i5k,a denotes the excitation energy of the quasi- 
particle in the a-th band and is given by -Ek.a — 

fi and Ak.a are the band 



?k,a + Ak^Q- Ck,Q — ek,a 

energy and the pairing gap of the a-th band, fi is the 
chemical potential. O^^j; = J2uy U^^.^yUu,^'/ is 
the matrix element of 1°' in the basis of the band eigen- 
state at momentum k. J7k,iy.a is the a-th eigenvector 
of the band Hamiltonian at momentum k. In deriving 
Eq.(l5]), we have used the inversion symmetry of the sys- 
tem and the fact that ZJ^ u' — v 

For comparison, we also calculate the Pauli spin sus- 
ceptibility of the band electrons. The Pauli susceptibil- 
ity is defined through the following spin-spin correlation 
function, 

X5(q,r) = -(g5MB)'(T,5^(q,r)5^(~q,0)), (6) 

where f/sA*B — 2 and the spin density operator at lattice 

site i can be written as = J2fi,s,s' '^i,^i,s'^s,s'^i,tJ.,s' ■ The 
uniform bare spin susceptibility can be shown to be given 
by 
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Ek^a — £'k+q,Q 

k+q,a ^ Ak,aAk-|-q,a 
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The difference between Eq.([5|) and Eq.Q lies in the 
fact that the orbital magnetic response has contribu- 
tion from both the intraband and interband processes, 
while the spin susceptibility has only intraband contri- 
bution, the reason for the latter is that the spin den- 
sity operator is diagonal in the orbital space. Such a 



difference will have important consequence in the su- 
perconducting state. As the coherence factor, (1 — 
Ck,QCk+q,c+ k,Q k+q,a \ vauishcs when o — > in the su- 

-fc'k.oi'k+q.a " ^ 

perconducting state, Xs{T) will show activation behav- 
ior if there is a full gap on the Fermi surface. On the 
other hand, the coherence factor for the orbital magnetic 

^k,aCk+q,,a+Ak.a^k+q./3 ' 



susceptibility, (1 



-), is in general 



-Ek.Q-Ek+q./3 

nonzero for a ^ /S when q — > 0. As a result, Xl(^) '^il^ 
in general be nonzero in the zero temperature limit. How- 
ever, the intraband contribution to the orbital magnetic 
susceptibility should be suppressed in exact the same 
manner as the spin susceptibility in the superconduct- 
ing state. 

Now we present the numerical results for both the or- 
bital and the spin magnetic susceptibilities. In our cal- 
culation, the chemical potential is self-consistently deter- 
mined by solving the particle number equation at each 
temperature. We then use Eq.® and Eq.Q to calculate 
the magnetic susceptibilities. The temperature depen- 
dence of the pairing gap is modeled by 



A(T) = AWl- 



T 



in which we have set A = 0.02eF and fc^Tc = 0.005eF. 
These parameters are typical for the Iron-based super- 
conductors. The band filling will be fixed at n = 6 at 
first, as the band structure calculation leading to model 
Hamiltonian Eq.([T]) is done at such a commensurate fill- 
ing. 
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FIG. 2: The spin and orbital magnetic susceptibility of 
the Iron-based superconductor calculated from the five-band 
model as functions of temperature. The band filling is fixed 
at n = 6. 

The results for the orbital and spin magnetic sus- 
ceptibilities are shown in Figl^l Both susceptibilities 
are found to show little temperature dependence below 
ksT = O.OSey in the normal state. This is in accordance 
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with the expectation for a typical band metaL The or- 
bital magnetic susceptibility is found to be several time 
larger than the Pauli spin susceptibility and is already 
comparable in magnitude with the observed total sus- 
ceptibility. More specifically, xf, {T) is found to be about 
3 times larger than the Pauli spin susceptibility xsiT), 
which is estimated to be about one-fourth of the observed 
total magnetic susceptibility at 150 K^^. The orbital 
magnetic susceptibility also show large anisotropy. The 
in-plane orbital susceptibility Xl (-^) found to be about 
1.8 times larger than the out-of-plane orbital susceptibil- 

ity xiiT). 

In our calculation, we have neglected the interaction 
correction. The interaction correction is believed to in- 
duce enhancement of the effective mass and thus the spin 
susceptibility. It is also believed that the interaction cor- 
rection will induce temperature dependence in spin sus- 
ceptibility. However, the interactions, such as the local 
Coulomb repulsion and the Hund's rule coupling, are not 
expected to renormalize the orbital magnetic susceptibil- 
ity directly. Thus, although the bare spin susceptibility 
calculated in this paper may not be a reliable estima- 
tion for the real spin response, the results for the orbital 
magnetic susceptibility should be robust. 

In the superconducting state, the spin susceptibility 
Xs {T) drops abruptly as the pairing gap develops on the 
Fermi surface and approaches zero in the zero tempera- 
ture limit. The in-plane orbital magnetic susceptibility 
x£(T) also exhibits a significant reduction in the super- 
conducting state, but remains nonzero in the zero tem- 
perature limit. On the other hand, the signature of the 
superconducting transition in Xl (^) almost unobserv- 
able. To understand why there is such a difference be- 
tween x£(T) and x%{T), we note that the diagonal ma- 
trix element O^. q, is identical zero as the result of the 
tetragonal symmetry of the system, while ^ ^ is in 
general nonzero. As we have mentioned above, the intra- 
band contribution to the orbital magnetic susceptibility 
also suffers from a suppression by the coherence factor in 
the superconducting state. 

To see how our results depend on the band filling, we 
have carried out the calculation for several different elec- 
tron concentrations. It is found that the qualitative fea- 
tures of both the spin and the orbital magnetic suscepti- 
bilities, for example, their temperature dependence and 
anisotropy, is quite robust, but their magnitudes are re- 
duced as we increase the electron concentration of the 
system. The detailed dependence of the susceptibilities 
and their relative ratio in the normal state on the band 
filling are shown in Fig. 3. The calculation is done at a 
fixed temperature fcsT = O.OSey and the rigid band ap- 
proximation is assumed. It is found that both the mag- 
nitude and the anisotropy of the orbital magnetic sus- 
ceptibility do not change significantly with band filling. 
On the other hand, the spin susceptibility show much 
stronger band filling dependence. 

Now we turn back to the experimental data shown in 
Fig. [T] As a common feature of all the data shown, both 
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FIG. 3: (a) The spin and orbital magnetic susceptibility as 
functions of the band filling at fixed temperature ksT = 
0.03el/. (b) The anisotropy ratio Xl/Xl as a function of 
band filling at fixed temperature ksT — O.OSeV . 



the in-plane and the out-of-plane susceptibility exhibit 
linear temperature dependence above Tc{or T^r) with al- 
most the same slope. The in-plane susceptibility is seen 
to be much larger than the out-of-plane susceptibility in 
all measurements. Thus the total susceptibility can be in- 
terpreted as consisting of two contributions: an isotropic 
component with a linear temperature dependence and a 
temperature independent component that is anisotropic. 
It is quite natural to attribute the isotropic component to 
the spin magnetic response, which with interaction cor- 
rection can exhibit strong temperature dependence22ri2i. 
The remaining anisotropic component is more likely to be 
the contribution of the orbital magnetic moment, which is 
not directly renormalized by the usual local electron cor- 
relation effect and should be temperature independent in 
the normal state. In principle, the anisotropy in magnetic 
susceptibility can also be caused by the spin-orbital cou- 
pling effect. However, in the Iron-based superconduc- 
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tors, the spin-orbital coupling is negligible small. There- 
fore, we feel our proposal for the anisotropy is more re- 
alistic for the Iron-based superconductors. 

In summary, we have shown that the orbital angu- 
lar momentum of the conduction electrons in the Iron- 
based superconductors can play a significant role in its 
low energy physics. It contributes a large temperature- 
independent anisotropic component to the magnetic sus- 
ceptibility. As a result, the total magnetic response in 
the normal state can be separated into a temperature 
dependent isotropic part XsiT) and a temperature in- 
dependent anisotropic part xl{0), where 6 is the angle 



between the magnetic field and the normal of the Fe-As 
plane. In other words, a fit of the experimental suscepti- 
bility to the formula x(r, 9) = Xs{T) + Xl{0) should be 
feasible. We note that the orbital magnetic moment can 
also contribute to the relaxation of the nuclear spins and 
other fluctuation effect at low energy. A full investiga- 
tion of the dynamical orbital response will be presented 
in future works. 
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